Study of the branching ratio of -0(3770) — > DD in e+e — >■ DD scattering 
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Based on the data of BES and Belle, the production of DD in the e + e~ — ¥ DD scattering process 
is studied in this paper. We analyze the continuum and resonant contributions in the energy region 
from 3.7 to 4.4 GeV. In the \ 2 fit to data, we obtain the resonance parameters of -0(3770), the 
branching ratio of i/>(3770) — > DD decay by confronting the data to the theoretical formula where 
both the contributions of the resonances, continuum and interference effects are included. We 
obtain the branching ratio of i/>(3770) — > DD decay is 97.2% ± 8.9%, as well as the branching ratio 
of V(4040), V(4160) -»■ DD decays. 

PACS numbers: 13.25.Gv, 13.66.Bc, 14.40.Gx 



Physics of e + e~ annihilation at the energy region of 
3 — 5 GeV is interesting. It attracts the focus of both 
experimental and theoretical studies. The physics at 
this energy region involves several well-established vec- 
tor resonances, J/ip, ip(2S), 0(3770), -0(4040), V(4160), 

which are bound states of quark and antiquark pair cc. 
Studying the production and decays of these resonances 
can deepen our knowledge about dynamics of interactions 
between quarks. The resonance ?/>(3770) has a mass just 
slightly above the threshold of DD pair production. Its 
decays evade the suppression of the Okubo-Zweig-Iizuka 
(OZI) rule [l|. This is consistent with the fact that the 
width of 0(3770) is about 2 orders larger than those of 
J/ip and ■0(25') j2|. The J/ip and -0(25) can only decay 
into non-DD final states, which is suppressed according 
to the OZI rule. It is believed that 0(3770) decays dom- 
inantly into DD pair. All the well-established non-DD 
decay modes of -0(3770) only show up with the branch- 
ing ratios at the order of 10 -3 — 10 -4 . The measurement 
of many other decay modes of -0(3770) only gives upper 
bounds, which shows that decay rates of these non-DD 
decay modes should be smaller than 10~ 3 or 10~ 4 [2, 3]. 
The sum of all these well-measured decay rates is at most 
at the order of several percent. 

In the e + e~ collider the properties of the resonance 
-0(3770) are measured through the scattering process 
e + e _ — > -0(3770) — > f, where / can be any final states 
like DD or any other hadrons. The branching ratio of 
0(3770) ->■ DD and -0(3770) ->non-DD can be derived 
from the measured scattering cross section of e + e _ — > 
DD and e+e~ -> hadrons. Both BES and CLEO-c Col- 
laborations measured the cross section of e + e~ — > DD 
at the center-of-mass energy E c .m. = 3773 MeV several 
years ago 0, HJ , and their results are in good agreement 
with each other. 

The CLEO-c Collaboration also measured the cross 
section of e + e~ — > 0;(377O) — > hadrons at -Ec.m. = 3773 
MeV 6ij. The difference between this and the cross 



section of e + e — » -0(3770) — > DD is found to be 
(-0.01 ± O.OSlo;^) nb, which_ indicates that the de- 
cay rate of -0(3770) to non-DD is tiny. However, the 
measurement of the BES Collaboration gives that the 
branching ratio of 7^(3770) -> DD is (85.5 ± 1.7 ± 5.8)% 
or (83.6 ± 7.3 ± 4.2)%, and the decay branching ra- 
tio of 0>(377O) to non-DD is (14.5 ± 1.7 ± 5.8)% or 
(16.4 ± 7.3 ± 4.2)% @, H, wnich is not consistent with 
CLEO-c's measurement. 

On one hand, a large branching ratio of -0(3770) to 
non-DD contradicts the fact that the sum of the branch- 
ing ratios of all the well-established exclusive non-DD 
decays is not large enough to give such a large decay rate 
for the inclusive decay mode. On the other hand, it is dif- 
ficult to understand such a large branching ratio for the 
non-DD decays theoretically. Calculations based on the 
method of QCD predict that the branching ratios of both 
exclusive and inclusive decays of -0(3770) to non-DD final 
states are very tiny, the sum of them is at most about 5% 
Therefore, the decays of 0(3770) -> DD ( or non- 
DD, an equivalent expression for one problem) is still an 
unsolved problem for both theory and experiment. 

At the e + e~ collider, the branching ratio of -0(3770) — > 
DD can be derived from the measured cross section of 
e + e~ — > DD. In this work, we reanalyze the exper- 
imental data of the e + e~ — > DD cross section in the 
center-of-mass energy region from 3.74 to 4.4 GeV mea- 
sured by the BES [lj| and the Belle [l3[ Collaborations. 
We will include not only the contribution of ^(3770) it- 
self, but also the contributions of the other resonances 
with masses near 3770 MeV. We also include the contin- 
uum contribution and its interference effect with the res- 
onances in the energy region 3.74 to 4.4 GeV. We finally 
find that the continuum contribution can be explained 
as an effect of the tail of -0(25), whose mass is about 40 
MeV lower than the threshold of DD pair production. 
By including the contributions of all the resonances near 
3770 MeV and the interference effects, the branching ra- 
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tio of -0(3770) — > DD we derived is apparently differ- 
ent from that of the BES experiment [7|, [8| . The contri- 
butions of the resonances with masses below and above 
that of 0(3770) are not included in analyzing the data of 
the e + e~ — > DD cross section by the BES Collaboration 
0, H| • The effect of these contributions is important for 
deriving the branching ratio of ?/>(3770) — > DD decay. 

In Ref. [IH, one of us analyzed the data of the e + e~ — > 
DD cross section measured by the BES and CLEO-c Col- 
laborations at £ c .m. = 3773 MeV 0,0). The formula 
describing the e + e~ — > D°D° or D + D~ cross section is 
derived 



D°D°,D+D-) = 
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where s = (pi + P2) 2 , Pi and P2 are the momenta of D 
and D mesons, respectively, ran is the mass of D or I? , 
D°, rrijp i the mass of the ith resonance, a — 1/137 is the 
electromagnetic fine-structure constant, and Q c — 2/3 is 
the electric charge of the c quark. The first term in the 
absolute- value squared F D jj(s) describes the continuum 
contribution. The second terms in the summation are 
the contributions of all the possible resonances, which 
are described by the Breit-Wigner form. The IYs are 
the total decay widths of resonances tpi's, and (pi's the 
relevant phases of the resonance contributions. g^^D 
is the coupling of the resonance ipi and DD, which is 
defined by 

(D(pi)5(p 2 )|V(p) i ) 
= -W^ddc W ■(Pi~P2)(2n) 4 6 4 ( P - Pl - P2 ), (2) 

where is the polarization vector of ipi, and A stands 
for the polarization state. f^ )i is the decay constant of 
the «th resonance ipi, which is defined by 



(3) 



Using the definition of the decay constant, the leptonic 
decay width of the resonance ipi is 



^QWf 2 



(4) 



With isospin symmetry, the coupling g^.j^Q and the 
continuum function F d q(s) are the same for both the 
production of D + D~ and D°D°. The difference of the 
cross sections of e + e~ to D°D° and to D + D~ is caused 
by the phase-space difference of D°D° and D + D~ . With 
the coupling g^^Di the branching ratio of the vector 
resonance ipi can be obtained 



(!) Bn(ipi -> D U D U , or D+D~) 
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(5) 

One can define the branching ratio of ipi — > DD as the 
sum of ipi — > D°D° and D + D~ , and then the branching 
ratio of ipi — > DD is 



BK{ipi -> DD) 
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D + D can be obtained from Eq.([TJ) 

a(e+ e - - DD) = * ( * = + ' '? = 

3 s 5 / 2 

te|%|^ e^f. (7) 
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With the expressions of the leptonic decay width and the 
branching ratios of ipi — > DD in Eqs. (jlj) and (|6l), one 
can reexpress the cross section of e+e~ — > DD in Eq.Q 
as 



a(e + e ^ DD) = ^ a 



^-^ q(m?i, ) Q^JT iT eeiBIZirn 5 ] 2 ., „ 
F DB (s) + ]T - 2 , ■ r e I ' ( g ) 



where the function g(x) is defined as and BIZi is the branching ratio of ipi — > DD. 

! The data of the cross sections of e + e~ — > mea- 

Sf(«) = I => ( 9 ) sured by BES [13] and Belle [13] are analyzed with the 

sj{x- 4m 2 30 ) 3 / 2 + (x- 4m 2 D± ) 3 / 2 formula in Eq.@. The BABAR Collaboration has 
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also made similar measurements before Belle's experi- 
ment [jjjj . and their data are consistent with Belle's, but 
for simplicity, we only use Belle's data in our analysis 
since they have higher statistics. The data require that 
the continuum term should be chosen as 



F m 2 



1>(3770) 



a 



(10) 



where Fq and a are parameters to be fitted. The value of 
the parameter a is found to be approximately the mass 
squared of 0(2S I ). Therefore the continuum term can 
be identified as the virtual contribution of the resonance 
ip(2S), whose mass is below the threshold of DD produc- 



tion. Then the continuum term is chosen to be 

Co 



V(2S) + * TO V-(2S) r V(2S) 



(11) 



where Co is the parameter to be fitted. 

The other resonances included in analyzing the 
e + e~ — > DD scattering cross section are 0(3770), 
G(3900), 0(4040), and -0(4160). Both the resonances are 
parametrized as Breit-Wigner form except for G(3900), 
for which the square root times a phase factor is used 
according to BABAR's finding [Hj]. The cross section of 
e + e _ — s- DD scattering in the energy range from 3.7 to 
4.4 GeV can be expressed as 



er(e + e 



DD) 
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g( m V,(4040)) 6 V r V>(4040)r ee 2g^2 "^(4040) 
a S — "1^,(4040) + (4040) (4040) 
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G(3900) 



3 m y,(4160) i<fe|2 



TO V(4160) + * TO i/H4160)r!/>(4160) 



(12) 



where Co, ci and the phases 4>iS are free parameters which 
will be fitted in the \ 2 n t to the experimental data. In 
our fitting the values of the following quantities are fixed 
and taken from the PDG Q: m v , (2 s) = 3686.09 ± 0.04 
MeV, r v , (2S) = 317±9 keV; the leptonic width of 0(3770) 
is r ee i = 0.265 ± 0.018 keV; to^ (4040) = 40 39 ± 1 MeV, 
r^,(4040) = 80 ± 10 MeV, the leptonic width of 0(4040) 
is r ee2 = 0.86 ± 0.07 keV; m^, (4160) = 4 1 53 ± 3 MeV, 
IV(4i60) = 103 ±8 MeV, the leptonic width of 0(4160) is 
r ee3 = 0.83±0.07keV. The quantities ra^^o), r^( 3 7 7 o), 
and the branching ratios of 0(3770) -» DD, 0(4040) -> 
DD and 0(4160) — » DD are set free and fitted from the 
data of BES [13 and Belle 0. The values of M G(3900) 
and ctg(3900) are varied and finally fixed, which can give 
the best fit to the data. 

Figure []] shows the nominal fit to the data of the cross 



section of e + e — > DD measured by BES [12| and Belle 
[l3j |. The data at 14 different energy points under the 
0(3770) peak between 3.73 and 3.80 GeV are from BES 
experiments by using the e + e~ scan [l2j ]. while the rest 
of the data at 27 different energy points above 3.80 GeV 
are from the exclusive initial state radiation (ISR) pro- 
duction of DD events from electron-positron annihilation 
at a center-of-mass energy of 10.58 GeV at Belle with an 
integrated luminosity of 673 fb _1 [13]. The data are cor- 
rected by the ISR HECl, 

and are all Born cross sec- 
tions. The solid curve is the best fit to the data by using 
the formula in Eq. (|12[) . The dashed curve is the contribu- 
tion of the resonance 0(3770), and the dotted curve is the 
continuum contribution from the tail of ip(2S). The very 
asymmetric line shape (dashed line in Fig. [I} is caused 
by the phase-space factor in Eq. flT " 



J 



The fit gives the following values for the parameters Cq, C\ and the phases: 

c = 8.75 ± 0.71 GeV 2 , c x = 1.00 ± 0.35 GeV 1/2 , 

= -2.63 ±0.09, cj>x = -1.89 ±0.33, 

02 = -2.14 ±0.14, 3 = 1.91 ±0.44, (13) 
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FIG. 1: The cross section of e + e~ — > DD. The dots with error bars are the data measured by BES P3] and squares with error 
bars are data from Belle The data are cross sections corrected by the ISR. The solid curve is the best fit to the data, the 
dashed curve is the contribution of the resonance 0j(377O), and the dotted curve is the continuum contribution of ip(2S). (a) is 
the cross sections in the whole energy range from 3.7 to 4.4 GeV, while (b) is the detail in the region of the resonance -0(3770). 



where the unit for 0Vs is radian. The results for the 
mass and width of 0(3770) and the branching ratios of 
V>(3770), 0(4040), and ^(4160) -> DD are 

»V(3770) = 37 76 ± 1 MeV, (14) 
I\/,(3770) = 28.5 ±2.1 MeV, (15) 

Sft(V>(3770) -> DD) = (97.2 ± 8.9)%, (16) 
F£(V>(4040) -> DD) = (25.3 ± 4.5)%, (17) 
BH(ip(4l60) -> DD) = (2.8 ± 1.8)%, (18) 

and the parameters for G(3900) are Mg(39oo) = 3900 ±20 
MeV and ctg(3900) — 52 ± 23 MeV which are fixed in the 
fit. In the above fit, there are totally 41 data points 
measured by BES and Belle and 11 free parameters are 
floated as shown in Eqs. (p~3|) ^ (fT8"| . The fitted quality 
is X 2 / n d = 1-06. 

In the fit, we find that the correlation between the 
phase and £ft(0;(377O) -> DD) is 0.73 which is 
large and mainly due to the component of the tail of 
ijj(2S). There is also a large correlation (~ —0.42) be- 
tween the parameter c\ for the structure G(3900) and 
the £ft(i/>(3770) -> DD). The structure of G(3900) had 
been sug gested in Ref. [l9[ and confirmed by BABAR 
data [15| . 

To test the significance of ip(2S) in the fit, a fit has 
been done without the contribution of ip(2S). We find 
the quality x 2 / n d = 1-5 (the \ 2 m this fit is worse by 14 
for 2 degrees of freedom by comparing to the x 2 value in 
the nominal fit) , while the value of x 2 / n d is 1.06 in the 
nominal fit. Furthermore, in the fit without the contri- 
bution of 0(25), we obtain the BTZ{ip(3770) -> DD) = 
(85.9 ± 7.0)%, which is significantly smaller than the re- 
sult obtained in the nominal fit. 



Figure [T] shows that the virtual contribution of -0(25) 
is very large (the dotted curve). It is even larger than 
the contribution of -0(3770) when the colliding energy is 
above the resonance region of -0(3770). The contribution 
of any resonance to e + e~ — > DD depends on both the 
coupling of this resonance with the virtual photon and 
with the DD pair. The coupling of the resonance with 
virtual photon can be described by its decay constant, 
which can be extracted from the measured leptonic de- 
cay width of the resonance. From the data of the leptonic 
width of V>(2S) andV>(3770) 0, one can obtain the decay 
constant of ip(2S) is f^ns) — 297 MeV, while the decay 
constant of -0(3770) is f^ 3770 ) = 100 MeV. This indi- 
cates that the coupling of -0(25*) with the virtual photon 
is approximately 3 times of that of -0(3770). One can also 
obtain the coupling of -0(3770) and ip(2S) with DD from 
the fitted result of B7£(<0(377O) -> DD) and the param- 
eter c with identifying c = 9^,(2S)DDQcU(2S) m ^(2S), 
which is indicated by the numerator of the second term 
of Eq.(|7|). The obtained couplings of the resonances with 
DD are 5</>(377o)dd = 12.8 and g 4 , {2 S)DD = 12.0, i.e., 
the couplings of 0(3770) and ip(2S) with DD are ap- 
proximately the same. Therefore the reason for the large 
contribution of ip(2S) to e + e _ — > DD scattering comes 
from the large coupling of -0(25*) with the virtual photon. 

In Fig. [TJ we show one of the best solutions. With the 
current data, multisolutions in the fit are possible since 
both the four phases and branching ratios are floated in 
the nominal fit. We investigate these effects, and find 8 
solutions with comparable fit quality However 6 of the 
8 solutions are with the £?7?.('0(377O) — > DD) less than 
70%, which are not consistent with the fact that no non- 
DD decay modes with significant branching fraction have 
been found in experiment. The sum of the branching 
ratios of all the well-established non-D£> decay modes 
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arc at most 2%-3% 0. Hereafter we discard the non- 
physical results. By assuming a constant width 1^(3770) 
for -0(3770) resonance in Eq. (fT2")l , we show two possi- 
ble physical solutions for comparison in Table |TJ While, 
by assuming an energy-dependent width for -0(3770) in 
the fit, we obtain similar results which are listed in the 
third column in Table [II The energy-dependent width 
r</,(3770) (s) is defined as 

r ^(3770)(s) = Td + D~(s) + ^D°D°( s ) + r non-flfi( s )>( 19 ) 

where T D+D - (s), T do qo{s) and T non _ D f)(s) are the par- 
tial widths for 0(3770) -> D+D~, 0(3770) -> D°D° and 
-0(3770) — mon-DD, respectively, which are taken in the 
form [7j 



D + D 



(s) = (3770)6* (E c 



2m D ±] 



(20) 



P D ± 



l + (rp° D± ) 2 
1 + (rp D ±) 2 



D°D° 



1 ^(3770) 

( Pd°\ 
{p%J 



9(E c . m . - 2m D o) (21) 
l + (rp° D o) 



1 + (rp D o) 



and 



^non-DDi s ) — ^ 



V>(3770) 



(1 - S+_ - S 



(22) 



where p 1 ^, and Pr> are the momentum of the D mesons 
produced at the peak of 0(3770) and at the center-of- 
mass energy a/s, respectively; r is the interaction radius 

of the cc, which is set to be 1.0 fm here; and Boo are 

the branching ratios for 0(3770) -> D+D~ and D°D°, 
respectively; and 6(E c . m . — 2m D ±) and Q(E c .m. — 2m£>o) 
are the step functions to account for the thresholds of DD 

production. In the fit, we fix the ratio Boq/B^ = 1.33, 

so that no additional free parameter is introduced. 

Our fitted result for the branching ratio -0(3770) — > 
DD is different from that of the BES Collaboration 
In our fitting we find that the contribution of the con- 
tinuum term from the tail of tp(2S), the contributions of 
the resonances 0(4040), 0(4160), the structure G(3900) 
and their interference effects with the resonance 0(3770) 
are important; they can seriously affect the fitting re- 
sult of the branching ratio of 0(3770) — > DD. If these 
effects are not included, a smaller result for the decay 
rate of 0(3770) — > DD will be obtained. The smaller 
decay rate of 0(3770) — > DD means a larger decay rate 
of 0(3770) — > non-DD. However, no large exclusive non- 
15 D decay mode of 0(3770) has been seen in experiment 
up to now. The sum of all the well-established exclu- 
sive non-DD decay rates is less than 2%-3%. This makes 
a puzzle for DD or non-DD decays of 0(3770). Our 
new analysis can solve this problem. With the branch- 
ing ratio of 0(3770) -> DD is (97.2 ± 8.9)%, the long- 
standing problem for 0(3770) decay disappears. Our re- 
sult can also be well understood theoretically. Theoret- 
ical prediction based on QCD to the branching ratio of 



Variables 


Constant width 


s-dependent 




Solution 1 


Solution 2 




™</>(3770) (MeV) 


3776± 1 


3776±1 


3780±1 


1^(3770) (MeV) 


28.5±2.1 


28.7±2.1 


29.7±1.3 


BKi (%) 


97.2± 8.9 


101.1±9.0 


98.3±10.4 


BTL 2 (%) 


25.3± 4.5 


34.7±4.8 


25.0±4.6 


B1Z 3 (%) 


2.8± 1.8 


40.4±3.8 


2.9±1.7 


CO 


8.75± 0.71 


8.67±0.67 


10.77±0.69 


Cl 


1.00±0.35 


0.82±0.29 


1.17±0.34 


0(rad.) 


-2.63± 0.09 


-2.56±0.09 


-2.49±0.08 


0i(rad.) 


-1.89± 0.33 


-1.55±0.36 


-2.32±0.30 


2 (rad.) 


-2.14± 0.14 


-1.62±0.11 


-2.56±0.21 


03(rad.) 


1.91± 0.44 


-3.03±0.1 


1.44±0.48 



TABLE I: Different solutions in the fit to data. The solutions 
1 and 2 are two physical solutions with constant width for 
V0377O) resonance. The last column is one solution for a fit 
to data with energy-dependent width (s dependent) as denned 
in Eq. (19]). 



the non-DD decay of 0(3770) is at most 5% [11]. This 
is consistent with our result of a large branching ratio of 
0(3770) DD decay [see Eq.flU])]. 

In summary, we analyze the data of the cross section 



of e + e 



DD measured by BES and Belle. The data 



are analyzed by including the contributions of all the res- 
onances in the energy region from 3.7 to 4.4 GeV. The 
contributions of the resonances and their interference ef- 
fects are important. Especially the virtual contribution 
of 0(2,5) is crucial for obtain our fitting results. Our re- 
sult for the branching ratio of 0(3770) — > DD can solve 
the problem long standing for 0(3770) — > DD and non- 
DD decays. We also get the branching ratios of 0(4040) 
and 0(4160) -> DD. 
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